NASA  Contractor  Report  18759i 
ICASE  Report  No.  91-51 


AD-A239  522 


ICASE 


STABILITY  OF  A  NON-ORTHOGONAL  STAGNATION 
FLOW  TO  THREE  DIMENSIONAL  DISTURBANCES 


D.  G.  Lasseigne 
T.  L.  Jackson 


Contract  No.  NAS  1-1 8605 
June  1991 

Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  Virginia  23665-5225 

Operated  by  the  Universities  Space  Research  Association 


DT1C 

ELECTS ‘ 

AUG  13 1991]] 


■u 


IUA5A 

National  Aeronautics  and 
Space  Administration 

Langley  Research  Center 

Hampton,  Virginia  23665-5225 


' "  PtaTRlBimON  St ATEMENT  A 

Approved  for  public  relecaei 
ljtotrfoutio®  Ualimitwd 


91-07634 


■ 


fjl  8  12  091 


STABILITY  OF  A  NON-ORTHOGONAL  STAGNATION 
FLOW  TO  THREE  DIMENSIONAL  DISTURBANCES 
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Abstract.  A  similarity  solution  for  a  low  Mach  number  non-orthogonal  flow  impinging  on  a  hot  or 
cold  plate  is  presented.  For  the  constant  density  case,  it  is  known  that  the  stagnation  point  shifts  in  the 
direction  of  the  incoming  flow  and  that  this  shift  increases  as  the  angle  of  attack  decreases.  When  the 
effects  of  density  variations  are  included,  a  critical  plate  temperature  exists;  above  this  temperature  the  stag¬ 
nation  point  shifts  away  from  the  incoming  stream  as  the  angle  is  decreased.  This  flow  field  is  believed  to 
have  applications  to  the  reattachment  zone  of  certain  separated  flows  or  to  a  lifting  body  at  a  high  angle  of 
attack.  Finally,  we  examine  the  stability  of  this  non-orthogonal  flow  to  self-similar,  three  dimensional  distur¬ 
bances.  Stability  characteristics  of  the  flow  are  given  as  a  function  of  the  parameters  of  this  study:  ratio  of 
the  plate  temperature  to  that  of  the  outer  potential  flow  and  angle  of  attack.  In  particular,  it  is  shown  that 
the  angle  of  attack  can  be  scaled  out  by  a  suitable  definition  of  an  equivalent  wavenumber  and  temporal 
growth  rate,  and  the  stability  problem  for  the  non-orthogonal  case  is  identical  to  the  stability  problem  for 
the  orthogonal  case.  By  use  of  this  scaling,  it  can  be  shown  that  decreasing  the  angle  of  attack  decreases 
the  wavenumber  and  the  magnitude  of  the  temporal  decay  rate,  thus  making  nonlinear  effects  important. 
For  small  wavenumbers  it  is  shown  that  cooling  the  plate  decreases  the  temporal  decay  of  the  least  stable 
mode,  while  heating  the  plate  has  the  opposite  effect.  For  moderate  to  large  wavenumbers,  density  varia¬ 
tions  have  little  effect  except  that  there  exists  a  range  of  cool  plate  temperatures  for  which  these  distur¬ 
bances  are  extremely  stable. 


This  work  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Contract  NASI -18605  while  the  authors 
were  in  residence  at  the  Institute  for  Computer  Applications  in  Science  and  Engineering,  NASA  Langley  Research  Center,  Hampton,  VA 
23665. 
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1.  Introduction.  A  similarity  solution  for  the  steady,  incompressible,  isothermal,  stagnation-point 
flow  impinging  on  a  flat  plate  at  an  angle  has  been  studied  independently  by  Stuart  (1959),  Tamada  (1979), 
and  Dorrepaal  (1986).  The  problem  is  solved  by  combining  an  irrotational  stagnation-point  flow  with  a 
linear  shear  flow  parallel  to  the  plate.  The  stream  function  splits  into  a  normal  and  tangential  component, 
and  it  is  discovered  that  the  normal  component  satisfies  the  same  non-linear  equations  which  govern  the 
orthogonal,  stagnation-point  flow.  The  tangential  component  of  the  stream  function  satisfies  a  linear  equa¬ 
tion  which  is  coupled  to  the  solution  for  the  normal  component.  The  major  consequence  of  adding  a 
tangential  component  to  the  flow  is  the  displacement  of  the  stagnation  point  toward  the  incoming  stream, 
which  increases  with  decreasing  angle.  Dorrepaal  also  discovered  that  the  ratio  of  the  slope  of  the  dividing 
streamline  far  from  the  plate  to  the  slope  of  the  dividing  streamline  at  the  plate  is  independent  of  the  angle 
of  attack. 


In  the  absence  of  a  shear  component,  Wilson  and  Gladwell  (1978)  have  shown  that  the  incompressi¬ 
ble,  planar  stagnation  flow  is  always  stable  to  three  dimensional  self-similar  disturbances.  Since  it  is  known 
from  experimental  observations  that  streamwise  vortices  may  be  formed  in  the  region  of  stagnation  flow 
past  bluff  bodies,  Lyell  and  Huerre  (1985)  re-examined  the  planar  stagnation  flow  problem  with  the  aid  of  a 
Galerkin  method.  In  particular,  they  verified  the  results  of  Wilson  and  Gladwell  and  characterized  higher 
stable  eigenvalue  branches,  showing  that  after  the  initial  branch  (found  by  Wilson  and  Gladwell)  the  other 
branches  come  in  pairs.  In  addition,  by  a  nonlinear  analysis,  Lyell  and  Huerre  showed  that  this  same  flow 
is  unstable  for  disturbances  of  sufficiently  high  amplitudes. 

Hall,  Malik,  and  Poll  (1984)  extended  the  work  of  Wilson  and  Gladwell  (1978)  to  include  both  the 
effect  of  crossflow  in  the  freestream  and  the  effect  of  suction  or  blowing  at  the  wall.  They  showed  that  the 
two-dimensional  flow  considered  by  Wilson  and  Gladwell  can  be  destabilized  with  sufficient  crossflow.  In 
particular,  with  no  suction,  the  three-dimensional  flow  becomes  linearly  unstable  at  a  critical  Reynolds 
number  (based  on  the  crossflow  velocity  in  the  freestream)  of  583.1.  The  effect  of  suction  and  blowing  is  to 
stabilize  or  destabilize  the  flow,  respectively.  The  weakly  nonlinear  analysis  and  full  numerical  solution  of 
the  three-dimensional  flow  is  given  in  Hall  and  Malik  (1986). 


It  is  the  purpose  of  this  paper  to  explore  the  effects  of  density  variations  on  a  two-dimensional  non- 
orthogonal,  stagnation  flow  and  then  to  study  its  stability  to  three  dimensional  disturbances.  In  Section  2,  a 
similarity  solution  is  obtained  assuming  small  Mach  number  and  an  ideal  gas.  A  stream  function  can  be 
defined  which  is  split  into  its  normal  and  tangential  components  and  is  required  to  satisfy  asymptotic  condi¬ 
tions  similar  to  the  conditions  proposed  by  Dorrepaal.  The  corresponding  differential  equations  governing 
the  temperature  of  the  fluid  and  the  normal  and  tangential  components  of  the  flow  field  arc  obtained.  In 
Section  3  the  formulation  of  the  stability  problem  is  given.  Section  4  contains  a  presentation  of  our  results 
of  the  stability  problem,  and  conclusions  arc  given  in  Section  5. 

2.  The  Mean  Flow.  It  is  well  known  that  the  continuity  equation  for  an  incompressible  fluid  can  be 
exactly  satisfied  by  introducing  a  stream  function 
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¥  =  P*;y  ,  ¥  =  y2 .  (2-2) 


which  represent  orthogonal  stagnation-point  flow  and  linear  shear  flow,  respectively.  The  overbar  represents 
a  dimensional  variable,  and  (5  is  the  dimensional  rate  of  strain.  The  solution  found  by  Stuart  (1959), 
Tamada  (1979),  and  Dorrepaal  (1986)  involved  any  linear  combination  of  these  two  flows.  The  combination 
chosen  here  is 

\g  =  ( y  y 2  cos  8  +  x  y  sin  8), 

(2.3) 

where  8  is  related  to  the  angle  of  attack  0A  by 

tan  8a  =  2  tan  9. 

(2.4) 

The  streamlines  given  by  the  family  of  curves 

y  =  -x  tan9  ±  Vjc2tan20  +  2¥p_1  see 0 

(2.5) 

are  shown  in  Figure  1  for  0  =  n/4.  The  dividing  streamline  is  found  by  setting  \|/  = 

m  =  -2  tan  0. 

0  and  has  slope 

Using  (2.1)  and  (2.3)  the  solution  to  the  potential  flow  is  now  given  by 

0  =  j$(y  cos0  +  it  sinQ),  V  =  -pysin0, 

(2.6) 

f  \ 

v2  v2 

P  =  P~ -p^^sin2©  Y  +  2  * 

(2.7) 

f  =  f .  , 

(2.8) 

where  P  is  the  dimensional  pressure  found  from  the  momentum  equations,  and  T  is  the  dimensional  tem¬ 
perature  assumed  to  be  everywhere  a  constant.  The  solution  for  the  constant  density,  but  varying  tempera¬ 
ture  profile,  has  recently  been  given  by  Lycll  (1990).  Although  the  streamfunction  (2.3)  satisfies  the 
Navier-Stokes  equations,  it  docs  not  satisfy  the  no-slip  conditions  at  the  plate  y  =  0.  For  this  reason  a  boun¬ 
dary  layer  must  exist  near  the  plate,  and  in  the  usual  fashion,  is  described  below.  The  above  potential  flow 
field  and  its  corresponding  boundary  layer  is  believed  to  have  applications  to  the  reattachment  zone  of  cer¬ 
tain  separated  flows  or  to  a  lifting  body  at  a  high  angle  of  attack. 

Following  the  analysis  of  Stuart  (1959),  Tamada  (1979)  and  Dorrepaal  (1986),  solutions  of  the 
Navier-Stokes  equations  in  the  limit  of  zero  Mach  number  are  sought  of  the  form 

U<*,r\)  =  xf'W>  +  4=  £'(fl),  (2-9) 


P  V  =-/0i). 


(2.10) 
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<,(">= m- 

f‘-T*fe-ikx  +  -£;Fw>-  (212) 

where  ti  is  given  by  the  Howarth-Dorodnilzyn  transformation 

y 

t|  =  J  p(y')dy'-  (2.13) 

0 

Here  the  temperature  T,  density  p,  and  pressure  P  have  been  non-dimensionalized  by  their  values  at 

infinity,  while  lengths  are  referred  to  a  characteristic  length  L  and  velocities  to  (5l.  In  addition,  y  and  V 

have  been  scaled  by  '/Re,  Re  being  the  Reynolds  number  of  the  flow.  The  viscosity  p.  is  assumed  to  be  a 
function  of  the  temperature.  For  small  Mach  number,  P  =  1  +  yM2P  so  that  P  represents  the  pressure 
deviation  from  the  ambient  The  Prandtl  number  Pr  is  taken  to  be  a  constant.  The  expansions  (2.9)-(2.12) 
are  required  to  match  with  the  outer  potential  flow  (2.6)  and  (2.7)  as  T|  °°  and  to  satisfy  the  no  slip  con¬ 
ditions  at  the  plate. 

Substituting  (2.9)-(2.12)  into  the  governing  equations  yields  the  following  system  of  equations  and 
appropriate  boundary  conditions  for  the  normal  component  of  the  flow  /,  the  temperature  T,  and  the 
tangential  component  of  the  flow  g : 

(j:f"Y  +  ff"  ~  (f')2  +  T  sin20  =  0,  / (0)  =  /' (0)  =  0,  /'(«)  =  sin  0,  (2.14) 

(^ry  +  Pr/r  =o,  r(0)  =  7>,  r(«)  =  i,  (2.15) 

(jg"Y  +  fg”  -r  8  +ATY~S m6cos0  =  O,  g  (0)  =  g'  (0)  =  0,  g’  (»)  =  q  cos  6.  (2.16) 

Here  TP  is  the  nondimensional  temperature  of  the  plate,  defined  as  the  ratio  of  the  dimensional  plate  tem¬ 
perature  to  that  of  the  outer  potential  flow,  and  A  =  A  cos  0  Vsin  0,  A  being  the  displacement  thickness 
whose  value  will  depend  upon  the  plate  temperature,  Prandtl  number,  and  the  particular  viscosity  law  under 
consideration.  When  7>  =  1,  the  above  equations  reduce  to  the  constant  density  case  with  A  =  0.6479r. 


f(r\)  =  a  F  (a  r\). 

(2.17a) 

T(  ti)  =  f  (aq). 

(2.17b) 

.  .  COS0  . 

gOl)-  ,  G(ati), 

(2.17c) 

The  change  of  variables 
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where  a  =  Vs  in  9,  reduces  the  above  set  of  equations  to  a  system  which  is  independent  of  the  angle.  The 
equations  in  their  simplest  form  are 

(4  F" )'  +  FF"  -  (F')2  +  f  =  0.  F  (0)  =  F'  (0)  =  0,  F' (<*>)  =  1,  (2.18) 

T 

(-4-  TO'  +  Fr  F  f'  =  0,  f(0)  =  7_,  f(oo)=l,  (2.19) 

T  p 

(ijrG"Y  +  FG"-FG'+Af  =  0.  G  (0)  =  G'  (0)  =  0,  G'(~)  =  q.  (2.20) 

Examination  of  the  above  equations  reveals  that  the  tangential  component  of  the  flow  does  not  affect 
the  normal  component  of  the  flow  nor  the  temperature  of  the  fluid.  In  fact,  equations  (2.18)  and  (2.19)  are 

identical  to  the  governing  equations  for  an  orthogonal  flow.  In  Figure  2,  F',  7  ,  and  G'  are  shown  for  the 

particular  case  of  the  linear  viscosity  law  p  =  f,  7>  =  0.8,  and  unit  Prandtl  number  (for  Chapman’s  linear 
law  p  =  C  T,  the  constant  C  can  be  scaled  out  by  rescaling  q,  F,  G,  and  A  appropriately,  but  must  be 
borne  in  mind  when  transforming  the  variables  back  to  their  dimensional  forms). 

A  quantity  of  importance  is  the  viscous  displacement  of  the  stagnation  point.  Note  that  as  q  -»  0, 

F  =  q2  -  ~2~ q3  +  •  •  •  ,  G  =  ^-q2- ^q3  +  •  •  •  ,  (2.21) 

2  6  2  3 

where  C  and  D  arc  chosen  numerically  so  as  to  satisfy  the  boundary  conditions  at  infinity.  For  the  constant 
density  case,  Tp  =  1,  C  =  1.232588,  and  D  =  1.406544.  Cx  and  £>,  are  known  constants  found  by  substi¬ 
tuting  the  expansions  into  the  original  equations.  The  viscous  displacement  of  the  stagnation  point  can  be 
found  by  writing  (2.21)  in  terms  of  /  and  g , 

/  =  ^-T|2(sin0)3/2  -  — -  rj3  sin2  8  +  •••,  g  ~  -^-q2cos0  — q3cos0  VsiinO  +  •••.  (2.22) 
2  6  2  3 

This  leads  to  the  following  expansions  of  the  velocity  field  near  the  plate 

U  ~  x  (C  1\ - n2  Vsin  9 )( sin 8  )3/2  +  (D  n  cos 9  -  D  j  rt2 cos 0  Vsin  0),  (2.23) 

2  VRe 

V  =  -  —  TP  ( sin  0  )3/2  q2 .  (2.24) 

The  dividing  streamline  near  the  plate  is  found  by  setting  (2.23)  to  zero;  then,  taking  the  limit  q  — »  0  pro¬ 
vides  the  following  expression  for  the  viscous  displacement  of  the  stagnation  point 

vkfcos0( sin 0 )~m.  (2.25) 

For  the  constant  density  case,  the  ratio  D  1C  =  1.141131.  The  stagnation  point  shifts  in  the  direction  of  the 
incoming  flow,  and  this  shift  increases  as  the  angle  decreases  from  re  /  2.  In  Figure  3  the  ratio  D  1C  versus 


Tp  is  shown  using  the  linear  viscosity  law  and  unit  Prandtl  number.  Note  that  a  critical  plate  temperature 
exists  around  7.4  for  which  the  ratio  vanishes.  For  plate  temperatures  above  this  value  the  stagnation  point 
shifts  away  from  the  incoming  flow  as  the  angle  is  decreased. 

3.  Formulation  of  the  Stability  Problem.  The  problem  considered  here  is  the  stability  of  the  non- 
orthogonal  stagnation  flow  described  in  the  previous  section.  Rescaling  the  transverse  coordinate  z  and  velo¬ 
city  W  in  the  same  manner  as  y  and  V ,  respectively,  and  keeping  only  the  leading  order  terms  in  a  large 
Reynolds  number  expansion,  it  is  easily  shown  that  the  non-dimensional  equations  which  govern  the  stabil¬ 


ity  of  this  flow  are  given  by 

P,  +  (PC/),  +  (pV)y  +  (p  W\  =  0,  (3.1a) 

p  [U,  +  UUX  +  VUy  +  WUt\  +  Px  =  biUy)y  +  Gi  Uz\ ,  (3.1b) 

p[V,  +  UVX  +  Wy  +WVt]  +  R  Py  =  (yiUy)x  +  (2\iVy)y  +  [p(V,  +  W,)]t  -  Ky,  (3.1c) 

p [W,  +  UWx+VWy  +  WWt)  +  RP,=  (p£/,),  +  (2 \iWX  +  [p(V,  +  Wy)]y  -  K, .  (3. Id) 

PIT,  +  UTX  +  VTy  +  WTI]  =  Pr~l  (pi,),  +  Pr~x  (pTJ, ,  (3.1e) 

pT  =  1,  p  =  7\  (3.1f,g) 


2 

where  K  =  —  p  (Ux  +  +  Wz ).  The  linear  viscosity  law  (3.1g)  has  been  used,  and  the  Mach  zero  limit  has 

been  taken.  Note  that  the  pressure  can  be  eliminated  from  (3.1c)  and  (3. Id)  in  favor  of  the  vorticity  com¬ 
ponent,  £2  =  -  Vz. 

Guided  by  the  fact  that  the  incoming  angle  can  be  scaled  out  of  the  mean  flow  equations,  solutions 
are  sought  in  the  following  form 

U (x,q,z  ,0  =  a2 x  F'(a  q)  +  a~:  Re"’/2cos  9G'(a  q)  +  e  [x  a2ut  (a  q)  +  Re~1/2«2(a  "H)]  EXP ,  (3.2a) 


V  (x  ,q  ,z  ,0  =  -  a  f  (a  q )F(a  q)  +  za  f(a  q)  v,(a  q )EXP,  (3.2b) 

VT(x,q,z,r)  =  ea  r\)EXP ,  (3.2c) 

T (x ,q,z ,i)  =  f(a  q)  (1  +  eT,(a  q )EXP),  (3.2d) 

p(x,q,z,r)  =  f~\a  q)  (1  -  eT,  (a  q)EXP),  (3.2e) 

£2(x,q,z,0  =  ea2f-\a  q)£2,(a  r\)EXP ,  (3.20 


where 
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EXP  =  exp{  i  a  a;  +  a2a>r  ),  (3.2g) 

and  e  <  1.  Note  that  the  size  of  the  pressure  disturbance  must  be  0(1/Re)  so  as  to  balance  equations  (3.1c) 
and  (3.  Id)  and  hence  does  not  effect  (3.1b).  Here,  a  is  the  span  wise  wavenumber,  taken  to  be  real,  and  to  is 
the  real  temporal  growth  rate.  The  assumption  of  to  being  real  is  consistent  with  previous  investigations  (eg., 
Wilson  and  Gladwell,  1978)  and  was  verified  numerically  by  allowing  to  to  be  complex.  In  every  case,  con¬ 
vergence  resulted  in  to  being  real.  Hall,  Malik  and  Poll  (1984)  obtained  complex  values,  but  only  with 
crossflow.  Presumably,  the  effect  of  crossflow  is  to  destroy  the  two-dimensional  symmetry  of  the  problem, 
thus  allowing  for  complex  eigenvalues.  The  disturbance  is  considered  amplified,  neutral,  or  damped 
depending  on  whether  co  is  positive,  zero,  or  negative,  respectively.  These  stability  results  are  limited  to  the 
special  class  of  self-similar  disturbances  considered  above.  Recently,  Brattkus  and  Davis  (1990)  considered  a 
wider  class  of  disturbances  for  the  stability  of  the  incompressible  Hiemenz  flow  and  concluded  that  these 
modes  were  more  stable  as  compared  to  the  self-similar  disturbances.  Therefore,  the  use  of  self-similar  dis¬ 
turbances  is  justified  for  the  constant  density  case,  and  there  is  yet  no  reason  to  believe  that  this  will  not 
also  be  true  for  the  variable  density  case,  although  this  still  needs  to  be  verified. 

Substituting  the  expansions  (3.2)  into  (3.1)  yields  the  linearized  equations  for  the  perturbation  ampli¬ 
tudes 

«,  +  v/  +  i  aw,  =  co 7’i  -  F  T\,  (3.3a) 

«i"  +  F  u{  -  (a2!2  +  co  +  IF' )«,  =  F"  v,  -  [2(F')2  -  IF  F”  -  7“]7',  -  F"  T (3.3b) 


fti  =  w\  -  i  af2v,,  (3.3c) 

0,"+  (F  -  //)£*,' -  [a 2f2+  0)  -  F'  +F  //  Ifij  =  d>,  (3.3d) 

~  T"  +  <F  +  4-ny 7Y  -  t  ~~  +  0)  +  2 F  //]  r,  =  //  v„  (3.3e) 

Pr  Pr  Pr 


where 


<t>  =  H  [  i  af2F  (1  +  2  Pr )  v  i  +  iaf2ui  -  ww,  j  -  2  i  a  f2F’  Tx’ 


-  i  at2 

V 


F"  -F  F’  +//[2u  +  4F'  -F2(l  +  2Pr)] 


T  i, 


(3.4) 


(3.5) 


Note  that  the  above  system  is  independent  of  the  angle  9  of  the  mean  flow,  and  hence  the  stability  problem 
for  the  orthogonal  case  is  identical  to  the  stability  problem  for  the  non-orthogonal  case.  The  above  system 
is  eighth-order  for  the  five  unknowns  uu  vj,  7j,  and  £2].  In  addition,  there  is  an  equation  for  u2  which 
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is  non-homogeneous  and  involves  G .  Once  the  system  (3.3)  has  been  solved,  u2  can  be  found  directly. 

The  appropriate  boundary  conditions  at  the  plate  are 

m1  =  Vi  =  h'1  =  7'i  =  0.  (3.6) 

Before  the  boundary  conditions  at  infinity  are  considered,  it  is  convenient  to  first  rewrite  (3.3)  as  a  real  sys¬ 
tem.  To  this  end  redefine  ux  =  Imag(.u  i),  vx=Imag{y  j),  wj  =  Real{w{),  T  ^  =  Imag(Tx),  and 
£l!  =  Real(£l{),  where  Real  and  I  mag  mean  the  real  and  imaginary  parts  of  the  complex  function  are  to  be 
taken.  Then  (3.3)  becomes 

“i  +  Vi '  +  ax'!  =  o)Tl  -  F  7Y,  (3.7a) 

II j"  +  F  II j'  -  (a2f2  +  CO  +  IF  )u,  =  F"  v,  -  [2 (F')2  -  IF  F”  -  f  ]T,  -  F"  T (3.7b) 


£2i  =  w,'  +  at 2  Vj, 


(3.7c) 


£2i "  +  (F  -//)«,'-  [  a2  f 2  +  co  -  F'  +  F//]Q,  =  O, 

+  (F  +  H)T\  -  [  +  to  +  2F  //]T,  =  //  v„ 

r  r  rr  rr 


where 


<t>  =  -// 


+  ar 


af2F  (1  +2 Fr) V[  +  a/2«i  +  cow,  j 
f  F"  -F  F’  +H  [2co  +  4F'  -F2(l  + 


+  2  a  7  2  F'  T\ 
2 Fr)]]  T,. 


The  asymptotic  behaviour  of  the  perturbation  amplitudes  at  infinity  are  given  by 

r1  =  X0(T)-A)<a2/'’'+<‘>)  +  e.s.f., 


(3.7d) 


(3.7e) 


(3.8) 


(3.9a) 


m  =X.,(ri-i4  )<“z  +  “  +  2>(i  +  <?(!))  +e.s.t.. 


(3.9b) 


where  e.s.t.  stands  for  exponentially  small  terms  as  r}  ->  «>.  Qj,  vIt  and  wj  have  similar  algebraic  terms  in 
their  asymptotic  expansions;  in  addition,  v,  and  w,  have  an  exponentially  growing  term  eaT1.  Wilson  and 
Gladwell  (1978)  have  shown  that  the  correct  boundary  conditions  for  the  stability  problem  is  that  all  pertur¬ 
bation  quantities  must  decay  exponentially  at  infinity.  Thus,  for  a  fixed  wavenumber  a,  the  eigenvalue  prob¬ 
lem  is  to  solve  the  eighth-order  system  (3.7)  subject  to  the  boundary  conditions  (3.6)  while  determining  a 
value  of  co  so  that  all  algebraic  and  exponentially  growing  terms  are  zero. 

Note  that  the  above  system  is  identical  to  that  considered  by  Wilson  and  Gladwell  (1978)  in  the  con¬ 
stant  density  limit.  They  solved  this  problem  by  use  of  invariant  embedding  and  Riccati  transformations. 
Lyell  and  Huerre  (1985)  solved  the  same  problem,  but  using  a  Galcrkin  method.  It  should  be  emphasized, 
however,  that  their  main  interest  was  the  nonlinear  stability.  Here,  the  above  eigenvalue  problem  is  solved 
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by  embedded  shooting  methods  using  a  Runge-Kutta  scheme  with  variable  step  size,  which  proved  to  be 
very  fast  and  very  inexpensive.  For  a  given  value  of  a,  an  initial  co  is  chosen  and  the  system  (3.7)  is 
integrated  subject  to  the  conditions  at  the  plate  (3.6)  and 

Q  i  =  Cj,  u\  =  C2 ,  n  =  C3,  Q',  =  C4.  (3.10) 

Then,  an  iteration  on  to  and  the  C,  ’s  by  a  Secant  method  is  performed  until  Xq,  A.i,  and  similar  constants  for 
the  other  quantities  were  less  than  10-6.  All  calculations  were  done  in  64-bit  precision,  taking  approximately 
one  second  of  CPU  time  on  the  CRAY-2  for  each  eigenvalue. 

4.  Results.  In  all  calculations  the  Prandtl  number  is  taken  to  be  one. 

As  discussed  in  the  Introduction,  for  the  constant  density  TP  =  1  case,  Wilson  and  Gladwell  (1978) 
have  shown  that  the  planar  stagnation  flow  is  always  stable  (co  <  0)  to  three  dimensional  self-similar  distur¬ 
bances.  In  their  paper  a  table  is  given  for  selected  values  of  the  wavenumber  a  with  corresponding  values  of 
the  growth  rate  co.  In  particular,  the  authors  showed  that  a  maximum  exists  in  to  for  an  a  of  about  0.298. 
In  addition,  they  suggested  that  there  were  other  more  stable  branches  than  were  shown,  but  this  was  not 
pursued  in  detail.  Lyell  and  Huerre  (1985)  re-examined  the  planar  stagnation  problem  with  the  aid  of  a 
Galerkin  method.  They  verified  the  results  of  Wilson  and  Gladwell  to  within  10%  using  a  one  term  Galer- 
kin  expansion  or  to  within  5%  using  a  two  term  expansion  and  displayed  two  more  eigenvalue  branches. 
The  authors  showed  that  after  the  initial  branch,  found  by  Wilson  and  Gladwell,  other  stable  branches  come 
in  pairs.  To  verify  the  present  numerical  scheme,  the  first  three  eigenvalue  branches  are  plotted  in  Figure  4. 
Results  for  the  first  branch  agree  with  those  of  Wilson  and  Gladwell  to  the  figures  they  gave.  The  other  two 
branches,  which  occur  as  a  pair,  agree  with  the  general  behaviour  given  by  Lyell  and  Huerre.  However, 
exact  agreement  is  not  possible  since  a  Galerkin  scheme  will  degenerate  for  higher  modes  unless  sufficient 
number  of  expansion  terms  are  kept.  The  calculations  presented  here  are  believed  to  be  correct  to  six  digits. 

The  purpose  of  this  section  is  to  extend  the  above  work  so  as  to  include  the  effect  of  variable  density. 
Results  for  the  least  stable  branch  are  plotted  in  Figure  5  which  shows  the  decay  rate  -  co  versus  the  plate 
temperature  TP  for  various  selected  values  of  a.  In  Figure  5a,  small  values  of  wavenumber  are  considered 
(a  <  0.4).  Note  that  for  a  =  0.001,  the  least  stable  branch  ceases  at  some  critical  plate  temperature  slightly 
greater  than  one.  As  a  increases,  so  does  this  critical  plate  temperature.  Most  importantly,  note  that  the 
effect  of  cooling  the  plate  is  to  reduce  the  value  of  the  least  stable  eigenvalue,  while  heating  the  plate  is 
seen  to  have  the  opposite  effect;  i.e.,  increasing  the  stability  of  these  modes.  If  the  plate  temperature  is 
increased  much  above  the  value  of  one,  then  each  of  these  small  wavenumber  modes  vanish  in  favor  of  a 
more  stable  branch.  The  plate  temperature  at  which  these  small  wavenumber  modes  vanish  increases  for 
increasing  wavenumbers.  In  Figure  5b,  moderate  to  large  wavenumbers  are  considered  (a  >  0.5).  For 
a  =  0.5  the  least  stable  branch  ceases  at  some  critical  plate  temperature  less  than  one,  about  0.3.  This  criti¬ 
cal  value  increases  as  a  increases.  Therefore,  it  is  seen  that  the  variation  in  density  has  little  effect  on 
moderate  to  large  wavenumber  disturbances  except  that  there  exists  a  range  of  cool  plate  temperatures  for 
which  these  disturbances  are  extremely  stable. 

It  is  now  clear  from  Figures  5a  and  5b  that  eigenvalues  on  the  least  stable  branch  can  cease  to  exist 
for  a  particular  combination  of  a  and  TP .  In  order  to  investigate  this  behaviour,  the  dispersion  relation  is 
shown  in  Figure  6  as  a  sequence  of  graphs  for  a  =  0.9  and  decreasing  values  of  the  plate  temperature. 
These  graphs  show  the  value  of  k,  versus  co  after  eliminating  all  other  algebraic  and  exponentially  growing 


terms.  Figure  6a  shows  the  dispersion  relation  for  the  case  TP  =  1.  First  note  that  for  to  >  1.5  the  curve  is 
increasing  and  approaches  zero  only  in  the  limit  co  -»  °°.  Thus,  never  vanishes,  and  the  flow  is  linearly 
stable.  For  to  <  1.5,  the  first  zero  of  X]  is  at  to  =  -1.853386,  which  agrees  with  Wilson  and  Gladwell  (1978). 
For  values  of  0)  less  than  this,  the  zeros  of  A.!  occur  in  pairs,  as  was  shown  by  Lyell  and  Huerre  (1985). 
Once  density  variations  are  taken  into  account,  no  matter  how  small  the  change  in  plate  temperature,  the 
structure  shown  in  Figure  6a  is  completely  altered.  Figure  6b  is  a  graph  for  TP  =  0.9.  Notice  that  the  verti¬ 
cal  asymptotes  have  changed  and  additional  ones  have  been  added  from  Figure  6a  to  6b.  Also  the  ordering 
scheme  of  Lyell  and  Huerre  for  the  constant  density  case  is  no  longer  valid.  Figure  6c  shows  that  some  of 
the  vertical  asymptotes  disappear  while  others  merge,  thereby  forming  a  continuous  curve  in  X,  over  a  broad 
range  of  to.  This  pattern  continues  through  Figures  6d  and  6e.  The  important  thing  to  note  is  that  between 
Figure  6d  and  6e,  the  least  stable  mode  has  disappeared. 

The  eigenfunctions  for  the  first  three  modes  of  Figure  6c  are  shown  in  Figures  7a,b,c,  respectively. 
The  first  two  modes  form  a  pair,  and  hence  their  eigenfunction  structures  are  similar.  The  third  mode  is  iso¬ 
lated,  and  its  eigenfunction  structure  is  different  from  those  of  the  first  two.  Also  note  that  each  eigenfunc¬ 
tion  decays  exponentially  at  infinity. 

The  effect  of  the  angle  0  of  the  mean  flow  on  the  linear  stability  can  now  be  assessed  from  (3.2g). 
For  any  given  value  of  TP ,  as  the  angle  is  decreased  so  is  the  wavenumber  and  the  magnitude  of  the  tem¬ 
poral  decay  rate,  thus  making  nonlinear  effects  more  important. 

4.  Conclusions.  It  has  been  found  that  variations  in  density  can  have  a  significant  feet  on  both  the 
structure  of  a  low  Mach  number,  non-orthogonal  stagnation  flow  and  its  linear  stability  to  three  dimensional 
self-similar  disturbances.  As  in  the  constant  density  case,  the  tangential  component  of  the  flow  can  be 
decoupled  from  the  normal  component  It  is  shown  that  the  viscous  displacement  of  the  stagnation  point 
due  to  the  non-orthogonal  flow  is  0{ ke~1'2),  which  depends  on  the  plate  temperature.  In  particular,  a  critical 
plate  temperature  exists  above  which  the  stagnation  point  shifts  away  from  the  incoming  flow  instead  of 
toward  the  flow  as  in  the  constant  density  case. 

For  the  linear  stability,  a  proper  scaling  can  be  defined  so  as  to  eliminate  the  angle  of  the  mean  flow 
when  considering  the  stability  characteristics.  Therefore,  the  stability  problem  for  the  orthogonal  problem  is 
identical  to  the  stability  problem  for  the  non-orthogonal  problem.  We  show  that,  for  small  wavenumbers, 
the  temporal  decay  rate  of  the  least  stable  mode  is  decreased  by  cooling  the  plate.  Finally,  we  show  that 
decreasing  the  angle  of  attack  decreases  the  magnitude  of  the  temporal  decay  rate,  thus  making  nonlinear 
effects  more  important. 

We  note  tha'  the  above  results  are  restricted  by  the  fact  that  the  similarity  solution  of  the  mean  flow  is 
not  appropriate  for  all  values  of  x.  We  are  currently  extending  the  above  analysis  to  a  fully  compressible 
subsonic  flow  by  means  of  the  IUingworth-Stewartson  transformation. 
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Figure  2a.  Plot  of  F'  versus  T)  for  TP  =  0.8  and  unit  Prandtl  number. 
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Figure  2b.  Plot  of  f  versus  t>  for  TP  =  0.8  and  unit  Prandtl  number. 
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Figure  6b.  Dispersion  relation  Xi  versus  w  for  a  =  0.9  and  TP  =  0.9. 
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Figure  6c.  Dispersion  relation  X[  versus  w  for  a  =  0.9  and  TP  =  0.8. 
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